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Abstract 

Using a narrow band positron beam, the response of a large high-resolution Nal(Tl) crystal to an incident 
positron beam was measured. It was found that nuclear interactions cause the appearance of additional 
peaks in the low energy tail of the deposited energy spectrum. 
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1. Motivation 

The PIENU experiment at TRIUMF Q is aim- 
ing at a measurement of the branching ratio R = 
r(7r ev -\- -K ^ ev^)lT{'K — >■ ^lv -|- tt ^ with 
precision <0.1%. The principal instrument used 
to measure positron energies from 7r+ e^u de- 
cays {E^+ = 70 MeV) and 7r+ — jjL^v followed by 
^+ e^uV decays (i?e+ =0 — 53 MeV) is a large 
single crystal Nal(Tl) detector Detailed knowl- 
edge of the crystal response is essential to reach- 
ing high precision, especially for determining the 
low energy tail response below 60 MeV 0]. In the 
following, results of measurements of the response 
of the Nal(Tl) crystal to mono-energetic positron 
beams are presented along with Monte Carlo (MC) 
simulations including photonuclear reactions. 



2. Experiment Setup 

The 48 cm diameter, 48 cm long Nal(Tl) crys- 
tal under study was surrounded by two adjacent 
rings of 97 pure Csl crystals Q. Each ring was 
comprised of two layers of 8.5 cm thick, 25 cm long 
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crystals. Positrons from the M13 beamline at TRI- 
UMF [i] were injected into the Nal(Tl) crystal to 
study its response. The positrons were produced 
by 500 MeV protons from the TRIUMF cyclotron 
striking a 1 cm thick beryllium target. After defin- 
ing the beam momentum at the first focus, the M13 
beam line is equipped with two more dipole mag- 
nets and two foci with slits before the final focus at 
the detector. The vacuum window was a 0.13 mm 
thick, 15 cm diameter Mylar foil. With this geome- 
try, slit scattering and the effect of the vacuum win- 
dow were expected to have negligible effect on the 
low energy tail. The incoming beam was measured 
with a telescope (see fig. [T|) consisting of 6 planes 
of wire chambers arranged in the orientation of X- 
U-V-X-U-V, where U(V) was at 60° (-60°) to the 
vertical direction, a plastic scintillator (5x5 cm^ 
area, 3.2 mm thickness), and the Nal(Tl) calorime- 
ter. The beam momentum width and horizontal 
(vertical) size and divergence were 1.5% in FWHM, 
2cm (1cm) and ±50mrad (±90mrad), respectively. 
The beam composition was 63% 7r+, 11% /i+ and 
26% e+. 



3. Measurement and Results 

A 70 MeV/c positron beam was injected into 
the center of the Nal(Tl) crystal. The beam 
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Figure 1: Schematic description of the experimental setup 
(not to scale) . The beam comes from the right and impinges 
on the Nal(Tl) crystal surrounded by two rings of 97 Csl 
crystals. In front of the Nal(Tl), there are 6 planes of wire 
chambers and a plastic scintillator. 

timing with respect to the 23 MHz cyclotron ra- 
dio frequency provided particle identification based 
on time-of-flight (TOF) together with the energy 
loss in the beam scintillator, allowing selection of 
positrons for studying the crystal response func- 
tion. Events due to positrons from decays of muons 
previously stopped in the Nal(Tl) crystal were sup- 
pressed by requiring wire chamber hits, and using 
TOF and pileup cuts. Pion and muon contamina- 
tion was reduced in the data to the 0.08% level. 

The Csl crystals were used in veto mode to se- 
lect events without shower leakage from the Nal(Tl) 
as well as for tagging events with delayed particle 
emission. Leakage from the NaI(Tl)'s downstream 
face was not detected but minimized by the 19 ra- 
diation length thickness of the crystal. 

The resulting positron energy spectrum is shown 
in fig. [5] (dark shaded histogram). The main peak at 
70 MeV has an asymmetrical shape due primarily to 
shower leakage with a width of 2.7% (FWHM). Sub- 
tracting the calculated beam momentum width in 
quadrature gave a Nal(Tl) crystal resolution of ap- 
proximately 2.2% (FWHM). Besides the main peak 
at 70 MeV, there are two additional structures at 
62 and 54 MeV. 

Studies were made to determine whether the ad- 
ditional peaks had either instrumental or physical 
origin. Using different settings of the momentum- 
defining and coUimating slits, which enhanced or 



suppressed slit scattering, no effect on the positron 
energy spectrum was found including the relative 
intensity of the peaks. Also, different tunes of the 
beamline {e.g. different focusing) did not change 
the measured energy spectrum. The beam momen- 
tum was varied in order to observe the correspond- 
ing position of the peaks. Fig. [2] also shows the 
spectra for 60 and 80 MeV/c beam momenta shifted 
and plotted on top of the reference histogram at the 
nominal momentum of 70 MeV/c. Signals from the 
Csl crystals were used to suppress the low energy 
tail due to shower leakage to enhance the second 
and third peaks. For all three beam momenta, the 
relative positions of the low energy peaks remained 
unchanged. The beam position dependence of the 
Nal(Tl) spectrum was also tested using wire cham- 
ber information, without finding any effect. Based 
on these tests, it is unlikely that there is an influ- 
ence of the beam settings in the appearance of the 
additional structures in the energy spectrum. 

In fig. [3] (top), the deposited energy in the 
Nal(Tl) crystal is shown as a function of the Csl hit 
time. The horizontal band at the beam energy cor- 
responds to accidental events, while the coincident 
ones from shower leakage are concentrated around 
ns. There are delayed events in the low energy re- 
gion that correspond to the second and third peaks. 
If delayed events between the vertical lines are se- 
lected, the shaded spectrum in fig. [3] (bottom) is ob- 
tained. The first peak (at approximately 70 MeV) 
is consistent with accidental coincidences. The sec- 
ond and the third peaks were enhanced after the 
delayed coincidence requirement. These results are 
consistent with the hypothesis of neutrons escap- 
ing the Nal(Tl) and giving a delayed signal in the 
Csl. Moreover, the energy deficits of the second and 
third peaks are consistent with the separation en- 
ergies for one {En — 9.14 MeV) and two neutrons 
{E2n = 16.3 MeV) emitted from ^'^''l. Since only the 
first hit is plotted in fig. [31 the observed secondary 
peaks are not due to the slow component of the Csl 
pulse. 

The yield of the second peak is consistent with 
the 30% solid angle and estimated 10% detection 
efficiency of the Csl calorimeter for neutron cap- 
ture. A delay of 100 ns is also consistent with the 
TOF of < 1 MeV neutrons. 

To estimate the number of neutrons involved in 
the second and third peaks, two Gaussian func- 
tions on a background with an exponential shape 
were fitted to both histograms in fig. [3] (bottom). 
The ratio N2 {N3) of the number of events in 
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Figure 2: Normalized Nal(Tl) energy spectra for incident 
positron beam momenta 60, 70, and 80 MeV/c. Tiie spectra 
were shifted and aligned to the peak at 70 MeV/c. His- 
tograms are scaled differently for easier comparison. 



the second (third) peak before and after the de- 
layed coincidence requirement is proportional to the 
product of the neutron detection efficiency and the 
number of neutrons involved, n2 (f^a)- The quan- 
tity = N3/N2 — n^/n2 indicates the ratio of 
the neutron multiplicities, which was found to be 

= 2.1 ± 0.2. Assuming that one neutron is in- 
volved in the second peak, this result suggests that 
the third peak arises when two neutrons escape 
from the crystal. 

The previous branching ratio experiment fsl was 
not able to detect these peaks because of the poorer 
energy resolution of the Nal(Tl) crystal (3 — 4% 
FWHM) employed. 

4. Simulation 

A MC simulation was developed, including all 
physics effects available in the GEANT4 package 
@i I3|- particular, photonuclear reactions with 
neutron(s) emission, scattering and absorption were 
taken into account using the QGSPJ3ERT physics 
processes list. In fig. SI the spectra obtained with 
the same detector setting with a monochromatic 
beam is shown. If in the simulation only electro- 
magnetic interactions were considered (dark shaded 
histogram), the low energy tail shows no structure. 
If hadronic interactions were included, additional 
structures appear (light shaded histogram), which 
are similar to those observed in the data (filled cir- 
cles). A closer look at the simulated data shows 
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Figure 3: (Top) Deposited energy versus Csl hit timing. 
(Bottom) The shaded histogram represents events selected 
by the timing cut (between the lines) shown on the top fig- 



that photonuclear reactions followed by neutron es- 
cape from the crystal are indeed responsible for the 
additional peak structures. Positrons entering the 
crystal produce an electromagnetic shower. One 
or more photons of the shower can be captured by 
^^^I nuclei. In the MC simulation, nuclear pho- 
toabsorption is generally followed by emission of 
neutrons (94%), protons (4%) or a-particles (2%). 
The kinetic energy and the separation energy of the 
neutron are not observed by the Nal(Tl) crystal if 
the neutron escapes. The second peak in the de- 
posited energy spectrum starts at i?i„ below the 
beam energy, where this reaction channel opens. 
According to the MC, the origin of the third peak 
in the spectrum is due to emission and escape of 
two neutrons. The neutrons can come from a sin- 
gle nucleus or from two separate ones (due to more 
photo-absorptions in the same shower). Both cases 
contribute to the third peak which starts at an en- 
ergy consistent with either the energy threshold of 
two neutron emission or twice the single separation 
energy Ein- 

The distribution of the kinetic energy for escap- 
ing neutrons is shown in fig. [S] The dashed his- 
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Figure 4: Comparison between data (filled circles with error 
bars) and simulation. The simulation was performed with 
(light shaded) and without (dark shaded) hadronic reaction 
contributions. The histograms are normalized to the same 
area. 

togram represents the kinetic energy of the neu- 
trons after nuclear emission, while the shaded his- 
togram shows the kinetic energy after escape from 
the Nal(Tl) crystal. The difference between the two 
spectra is due to elastic and inelastic scattering re- 
actions in the Nal(Tl) crystal. 

In fig. ini the correlation between the number 
of elastic scatterings and the neutron kinetic en- 
ergy at production is shown. Figs. [5] and [6] suggest 
that, although the primary source of the second and 
third peaks is low energy neutron emission from 
photonuclear reactions, many neutron elastic scat- 
terings significantly lower the escaping neutron ki- 
netic energy, returning "lost" energy to the Nal(Tl) 
crystal. 

The agreement between simulation and experi- 
ment is not perfect. Given the high number of in- 
teractions which a neutron can experience in a large 
crystal, a small error in the models for elastic and 
inelastic scattering can be amplified. Moreover, in 
GEANT4 photonuclear reactions are parameterized 
on a limited data set of nuclides. 

5. Conclusions 

The response of a large Nal(Tl) crystal to a 
positron beam of 70 MeV/c was investigated in 
preparation for the PIENU experiment at TRI- 
UMF. Low energy structures were observed in the 
energy spectrum and the mechanism for their origin 
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Figure 5: Simulation of the kinetic energy of the neutrons 
produced in (white histogram) and those that escaped from 
(shaded histogram) the Nal(Tl) crystal. 
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Figure 6: Simulation of the number of elastic scatterings as a 
function of the kinetic energy of the neutrons after escaping 
the nucleus. 

was found to be consistent with neutron emission 
due to photo-absorption followed by neutron escape 
from the crystal. 
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